Vibrational excitations in systems with correlated disorder 



W. Schirmacher and C. Tomaras 
Phys.-Dept. E13, Technische Universitdt Miinchen, Germany 

B. Schmid 

Fachbereich Physik, Universitdt Mainz, Germany 



o 
o 

(N 
> 

o 

oo 



C/3 



o 



> 
m 



o 



X 



G. Viliani 

Dipt, di Fisica, Universitd di Trento, Italy 
G. Baldi 

INFM-CNR CRS-SOFT OGG, c/o E.S.R.F., BeamUne ID16, BP220, F-38043 Grenoble Gedex, France 

G. Ruocco, T. Scopigno 
ipt. di Fisica, Universitd di Roma, Italy; CRS SOFT-INFM-CNR c/o Universitd di Roma, Italy 

We investigate a d-dimensional model (d — 2,3) for sound waves in a disordered environment, in 
which the local fluctuations of the elastic modulus are spatially correlated with a certain correlation 
length. The model is solved analytically by means of a field-theoretical effective-medium theory 
(self-consistent Born approximation) and numerically on a square lattice. As in the uncorrelated 
case the theory predicts an enhancement of the density of states over Debye's c^'*"^ law ("b oson 
peak") as a result of disorder. This anomay becomes reinforced for increasing correlation length 
^. The theory predicts that ^ times the width of the Brillouin line should be a universal function 
of ^ times the wavenumber. Such a scaling is found in the 2d simulation data, so that they can 
be represented in a universal plot. In the low-wavenumber regime, where the lattice structure is 
irrelevant there is excellent agreement between the simulation at small disorder. At larger disorder 
the continuum theory deviates from the lattice simulation data. It is argued that this is due to an 
instability of the model with stronger disorder. 

PACS numbers: 65.60 



I. INTRODUCTION 

The influence of quenched disorder on the dynamic 
properties of solids is enormous and is subject to 
widespread experimental and theoretical investigations 
The disorder leads to strong modifications of phys- 
ical properties of the solid. On the other hand, the ab- 
sence of lattice order in the solids and the correspond- 
ing breakdown of the Bloch theorem leads to appreciable 
difficulties for the theoretical interpretation of the dis- 
order-induced phenomena. Here mean-field theories and 
particularly effective-medium theories have been of much 
help [1], as they very often lead, at least, to a qualitative 
understanding of the influence of the disorder. In par- 
ticular the coherent-potential approximation [s^ (CPA) 
and its small-disorder version, the self-consistent Born 
approximation (SCBA) Q], have proved to be useful for 
interpreting the electronic and other spectral properties 
of disordered solids. In many cases physically different 
situations can be mapped onto each other. So one can 
convert an electronical problem to a vibrational one by 
replacing the energy E by — w^, where uj is the frequency 
parameter. If E is replaced by id;, one studies the math- 
ematical analogous diffusion problem 0, [^, Q , where uj 
denotes the time Fourier parameter of the diffusion dy- 
namics. 

In the case of vibrational properties the disorder-in- 



duced excess in the density of states (DOS) over Debye's 
u}'^~^ law {d is the dimensionality) ("boson peak") has 
been successfully explained for a lattice model by com- 
paring a simulation of a disordered lattice system with 
the predictions of the lattice CPA 0]. The comparison 
showed once more that the CPA is a reliable theory of 
disorder. In this study it was shown that the boson peak 
anomaly marks the crossover from plane-wave like vibra- 
tional states to disorder-dominated states with increasing 
frequency uj. Near the crossover the effective sound ve- 
locity V becomes complex and frequency dependent. This 
corresponds to a dc-ac crossover in the analogous diffu- 
sion problem Q . 

A similar but lattice-independent approach is the gen- 
eralised elasticity theory in which the disorder is assumed 
to lead to spatial fluctuations of the elastic constants 

8]. This model was solved by functional- integral tech- 
niques in which the SCBA plays the role of a saddle-point 
in a non- linear sigma- model treatment. This theory al- 
lowed for a generalization to include transverse degrees 
of freedom and to formulate a thermal transport theory 

9]. Within this theory it was shown that the so-called 
plateau in the temperature dependence of the thermal 
conductivity [lo| is caused by the boson-peak anomaly 
0. Furthermore, it has turned out [H] that the sound at- 
tenuation parameter (which is proportional to the imag- 
inary part oi v{uj)) is related to the excess in the DOS in 
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the boson peak regime. 

It should be noted that the model of spatially fluctuat- 
ing elastic constants, treated in CPA and SCBA is by no 
means the only approach for explaining the boson-peak 
anomaly. In fact, an enormous number of possible ex- 
planations have been published in the literature, which 
can roughly grouped into three classes: i) defect models, 
a) models associated with the glass transition and in) 
models with spatially fluctuating clastic constants. 

i): Defects with a heavy mass can produce resonant 
quasi-local resonant states within the DOS [l^, [H, [3 
and be thus the reason for the boson peak and the re- 
duction of the thermal conductivity. Similarly defects 
with very small elastic constants, near which anharmonic 
interactions are important (soft potentials), can produce 
quasi-local states, which, if hybridized with acoustic exci- 
tations may produce a boson peak [Tp . [TgI and a plateau 
in the thermal conductivity (TtI . llSj Inhomongeneities 
may also be the source of local vibrational excitations 
that contribute to the excess DOS Specifically 
in network glasses bond-angle distortions may also con- 
tribute to the boson-peak anomaly [13, [2l| . In a recent 
study the predictions of a defect model has been com- 
pared with those of a fiuctuating elastic constant model. 
a): In theories of the glass transition [H, [Hill [H [111 
the boson peak arises as a benchmark of the frozen glassy 
state. Hi) In models wit h q uenched disorder of elastic 
constants [1, [E i, i, [H [13 [H, [H, [M [IH the boson 
peak marks the lower frequency bound of a band of ir- 
regular delocalized states with random mutual hybridiza- 
tion. These states are neither propagating nor localized 
The models have been solved with the help of nu- 
merical simulations as well as effective-medium theories. 

A drawback of these models is that they are based on 
the model assumption of uncorrclatcd disorder, i.e. spa- 
tial fluctuations of the physical quantities are assumed 
to be uncorrelated. This assumption is only justified if 
the spatial correlation length is smaller or of the order 
of the natural length which appears in the system un- 
der consideration. In the present problem, namely vibra- 
tions in disordered solids, there are two important length 
scales. One is just the interatomic distance a, the other is 
the sound velocity, divided by the boson-peak frequency, 
which is in experimental data of the order of several in- 
teratomic distances. The wavenumber corresponding to 
this length scale is the maximum wavenumber wich can 
serve as a label for wave-like vibrational states. In any 
case it is a more sound procedure to start with a the- 
ory with correlated disorder and make the approxima- 
tion of short correlations (if appropriate) only in the end. 
Such theories are available [sH, [sj, [sj] . The aim of the 
present contribution is twofold: First we summarise the 
main features of the long-range-order SCBA. Secondly 
we present results of a two-dimensional simulation and 
compare them with those of the analytic theory. 



II. MODEL AND SELF-CONSISTENT BORN 
APPROXIMATION (SCBA) 

We start with the equation of motion for scalar wave- 
like excitations u(r,t) in a c?-dimensional disordered 
medium 



^2 

— M(r,t) - V - v^ir) • Vu(r,i) 



(1) 



Here v{r) is a sound velocity (and w^(r) an elastic con- 
stant), which is supposed to exhibit random spatial fluc- 
tuations v'^{r) — Wq + A(r) with Vq = (w^) and 



(A(r + ro)A(ro)) = C(r) 
C(k) - (A2)C7o[fc' + e 



^ c{k)^{A')Co[k'+r'r^ (2) 

with Co = 27r(3c?— 5)/^. In an effective-medium approx- 
imation the disordered system is mapped onto a homo- 
geneous system, in which the influence of the disorder 
enters via a self-energy function S(k, z) with z = uj + ie. 
The dynamic susceptibility is given by 

= -z2 + P(„2_s(k,z)) 

y function obeys the self- 



In SCBA the self-enere 
consistent equation [1, [sj 



S(q,z) = i J 



dk\ 

-) C(q-k)x(k,z) (4) 



|k|<fci 



with the Debye wavenumber kjj ~ [Zdir'^^^]^/^ /a. 

In the present study we are mainly interested in the 
low-frequency and -wavenumber properties, so we replace 
the self-energy by its g ^ limit S(z) and obtain the 
SCBA equation 



S(z) 



dkV 



k^C{k)/{A^) 



2 V>°^ J ^271 J -Z2 + fc2 (i;2 - I](2)) 
|k|<A:D 



(5) 



with the normalization 

/ (dk/27r)''C(k)/{A^> and the 

k|<fc_D 

The DOS is given by 



»H^^,J/(S 



constant tpd = 
"disorder parameter" 



^z^ + k^{v^-j:{z)) 



(6) 



In Fig. 1 we have plotted the "reduced DOS" g{uj)/Lu'^ 
for d = 3 and for three values of 7 and two values of 
£. First we notice that, as in the uncorrelated case [3, 
[1, [^ there exists a critical amount of disorder 7c — 0.5, 
beyond which the system becomes unstable. The "boson 
peak" becomes more pronounced and situated at lower 
frequencies as this value is approached. The interesting 
new feature in the correlated case is that the boson peak 
is re-inforced by the correlation and again shifted towards 
lower frequencies. 
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FIG. 1: Reduced DOS g{uj)/LO^ for ^ = l/kD (dashed lines) 
and ^ — b/ko (full lines) and for three disorder parameters 
(from left to right) 7 = 0.49, 0.47 and 0.45. 



FIG. 2: Scaled Brillouin linewidth F^/uo against scaled 
wavenumber for ^ = l/fco (dashed lines) [3^], C = ^/ko 
(full lines), ^—lO/ko and IS/fcu (dotted lines). The disorder 
parameters 7 are the same as in Fig. 1. 



III. DYNAMIC STRUCTURE FACTOR AND 
DENSITY OF STATES 

We now divide the self-energy function into a real and 
imaginary part E(2;) = S'(w ~ 0)-|-zS"(w) and define the 
renormalised sound velocity as v'^ = Vq—Ti'. The dynam- 
ical structure factor S'(k, lj), which can be measured by 
inelastic neutron or X-ray scattering, and which is the 
Fourier transform of the dynamic density-density correla- 
tion function, is then given by the fluctuation-dissipation 
theorem [33 as 



S{k,ij) 



knT 



TTLU 



Im{x(fc,z)} 



k^Y."{uj)/2uo 



27rw2 [{kv-ujy + (/s2I]"(tj)/2cj)2 



(7) 



The Brillouin resonance is given hy lj — kv and the line 
width (FWHM, sound attenuation parameter) is given 
by 



T{k) = -Y."{uo ^vk) 



(8) 

so that 



It can easily be shown that S" oc w'' for uo 
r oc 9'^"'"^ for g ^ (Rayleigh law). 

We now introduce the dimensionless variables 
t = Y,/vl, Co = CoC = 27r(3d- 5), i = z£,/vo=Cj + K 
and k = k^. We then obtain 



|k|<«fc 




d 



e [1 + 



-z2-f fc2(l_E(|)) 



(9) 



In this expression the correlation length ^ enters only via 
the upper k cutoff. For the limit ^ ^ a we therefore 
obtain the scaling relation 



f = re/«o = /(7,fc) 



(10) 



In Fig. 2 we have plotted F against A; for d = 3, for the 
three 7 values of Fig. 1 and for ^ = l/fcc, S/fc^, lO/fc/j, 
and 15/fc£). We see that the scaling is obeyed except for 
i = l/ko nil as expected. As in the uncorrelated case 
the boson peak (see Fig. 1) marks the crossover 
from the Rayleigh regime F cx A;^ to a behaviour F oc fc^ 
with s w 2. 



IV. SIMULATION 

We now discretize ([Ij in d = 2 on a square lattice, 
which then takes the form of an equation of motion for 
unit masses connected by springs with spring constants 



1 

2a2 



{v^{ri) + v^{rj)): 



t) 



E 

I n.N. 



(11) 



In the simulation [37| the force constants Kij are ex- 
tracted from a random distribution with mean Kq and 
variance cr^. The correlation is established following the 
Fourier filtering method (FFM) 3]. The network of 
random springs is created starting from a random set 
of {2N)'^ numbers uniformly distributed around zero ob- 
tained from a pseudo random number generator. The 
FFM method is then used to generate a 2N x 2N 
two-dimensional lattice of "pair" random numbers {fyy} 
which obey a spatial correlation 



-§(»+j)/« 



(12) 



The random spring constants with mean Kq and variance 
are extracted from the pair numbers rjij. The lattice 
spacing of masses-and-springs system is twice of that of 
the random-number lattice. Using this statistics the dy- 
namic structure factor S{k,uj) of the model has been de- 
termined by the method of moments [1^, SO, |4l| . In Figs. 
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FIG. 3: Symbols: Linewidths F of the simulated Brillouin spectra, multiplied with the correlation length ^ and divided by the 
sound velocity v. o : ^ = 2.79a, A : ^ = 5.42a, V : ^ = 10.56a, + : C = 20.55a, x : = 40.0a, Straight line: theory. The resuhs 
of the left picture correspond to a disorder parameter 7 = 0.04 that of the right picture to 7 = 0.111. 



3a and 3b the scaled widths of the Brillouin peak of sam- 
ples with different correlation lengths have been plotted 
against the scaled wavenumber k. It is clearly seen that 
the simulated data follow the predicted scaling law. For 
the case with the lesser disorder (7 = {cr/Ko)'^ ~ 0.04) 
there is very good agreement with the theory ( ^ with 
no cutoff in the integral) in the small wavenumber limit, 
where the lattice and continuum models should agree. In 
the high-wavenumber regime, of course, the lattice char- 
acter of the simulated system becomes distinct. 

Let us turn to the discussion of the data of Fig. 3b with 
the increased disorder 7 = 0.111. The continuum theory 
in this case predicts the Rayleigh law F oc fc'^ (continuous 
line). We checked the stability of the system by investi- 
gating the simulated density of levels g{u!'^) — g{ijj)l2uo 
and found that this quantity exhibits nonzero values for 
up < 0, which means that the system is unstable. For 
this case it is known that the imaginary part of the self 
energy E" is constant and passes continuously from pos- 
itive to negative values of in this case. Consequently 
the line width F(cj) oc kT," shows a linear increase for 
small omega. Such a behavior is obviously an artefact 
of constructing a harmonic model with too much dis- 
order, which leads to a small fraction of negative elastic 
constants. In a realistic physical system such "would-be" 



negative elastic constants are removed by the anharmonic 
interaction which causes relaxation of the system towards 
a stable situation. This has been nicely demonstratedby 
a model calculation of Gurevich et al flal. 



V. CONCLUSION 

We have investigated the vibrational properties of 
disordered systems with correlated disorder both analyt- 
ically by the self-consistent Born approximation as well 
as by a simulation applying the method of moments and 
the Fourier-filtering method. The sound attenuation con- 
stant F(fc) (width of the Brillouin line) is found to scale 
with the correlation length ^ in such a way that ^F is a 
universal function of (^k. The enhancement of the density 
of states (boson peak) is found to be re-inforced by the 
correlations. 
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